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Abstract 

We consider thermal QCD in the large Nc limit, mainly in 1+1 
dimensions. The gauge coupling is only taken into account to minimal 
order, by projection onto colour singlets. An expression for the free 
energy, exact as Nc — > do, is then obtained. A third order phase 
transition will occur. The critical temperature depends on the ratio 
Nc/L, where L is the (infinite) spatial length. In the high temperature 
limit, the free energy will approach the same value as in the free theory, 
whereas we have a mesonic like phase at low temperature. Expressions 
for the quark condensate, {^^}, are also obtained. 
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Considerable interest has recently been devoted to the assumed decon- 
finement phase transition in QCD at high temperature and/or density. Exact 
analytical results are incomprehensible sofar, and we have to rely on approxi- 
mate methods or computer simulations on the lattice. It has been well known 
for a long time that some insight in the confinement mechanism may be 
gained by considering the large Nc limit, for an SU(Nc) gauge theory. This 
limit is particularly fruitful in 1+1 dimensions (see e.g. Refs. |], @, 0, |], §])• 
QCD in 1+1 dimensions has also received attention recently due to the re- 
markable fact that different regularizations seems to lead to different models, 
even within the same choice of gauge. In his pioneering work, 't Hooft [1[] em- 
ployed a cutoff around the origin in momentum space, that later was shown 
to be equivalent to a principle value prescription . Within this framework, 
no free quarks can propagate, but they are confined into mesons, consisting 
of quark anti-quark pairs. Wu H], on the other hand, suggested another 
regularization, allowing for a Wick rotation to Euclidean space. The bound 
state equation is more complicated in this case, and no solution has yet been 
found. Bassetto and Nardelli with co-workers 0, 0, [| have compared the 
two regularizations, and computed for example the expectation value of the 
Wilson loop. The expectation value of the Wilson loop with Wu's regular- 
ization, was recently calculated by Staudacher and Krauth [[Hj , who showed 



that confinement is not enforced in this case, unlike using the 't Hooft reg- 
ularization. It has been suggested by Chibisov and Zhitnitsky (Tl| that the 
models due to different regularizations could show up as different phases, 
that makes a study of the system at finite temperature intriguing. In 2 di- 
mensional QCD, confinement is obvious as an infinite energy for coloured 
states. 

We shall in this letter take the interaction into account only by a projec- 
tion onto colour singlets, i.e. we neglect the coupling unless it multiplies the 
infra red divergence. The quark contribution to the free energy in 3+1 di- 
mensions is then easily obtained from the free energy presented here in 1+1 
dimensions through the substitution, L J dp/(2n) — > 2V j d 3 p/(2n) 3 . The 
projected free energy in 3+1 dimensions has earlier been considered in the 



high temperature limit by Skagerstam [12]. The full interacting case in 2 di- 
mensions, that is troubled by infra red divergences, will be more extensively 
considered in another article. Including the interaction in the strict large Nc 



limit L/Nc — > 0, McLerran and Sen |13|], and employing a different method 



Hansson and Zahed showed that in the low temperature 't Hooft me- 
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son phase, the partition function InZ = O(L), due to confinement. Using a 
different resummation scheme for L/Nq — > oo , McLerran and Sen |L3) also 



found another phase of almost free quarks, with In Z — O(LNc) ■ The phase 
transition itself was not obtainable, and could only take place as T — ► oo, 
i.e. T oc Nc- By considering finite « = Nc/L, we are here able to monitor 
the phase transition exactly, and for finite k, it will take place at finite T. 
However, due to the negligence of the interaction we do not have confinement 
in the strict sense. Our low temperature phase consists of quark — anti-quark 
pairs, showing up for example in Bose-Einstein distribution functions, but 
these are not the confined mesons of 't Hooft's, since InZ = O(LNc)- 

Let us now consider the unnormalized density matrix p = e~@ of free 
coloured fermions having the Hamiltonian 

H= f ^( P )h{p)^{p), (1) 

where 

h(p)= 1 °( ll p + m), (2) 

and Jp = J / 2 f„ , in n spatial dimensions. We shall here mainly discuss n — 1. 
In the functional representation the expression for p is jl5| 



PiViVh V*2V2) = N p exp((r]* + r^fy^i + r] 2 ) + V*iV2 ~ V2V1), (3) 

where 

N p = exp(AN c V /ln(cosh^o;)), (4) 
Jp 2 

(3 h 

Vt p = -tanh(-a;)-lc, (5) 
2 uj 



Lo{p) = \J P 2 + m 2 , (6) 

and L denotes the spatial length. In the expression for piViViiVZVz)] sum- 
mation over spinor, colour and momentum indices is understood. Moreover 
Nc denotes the number of colours and 1q is the unit matrix in colour space. 
We now wish to calculate the partition function for colourless states. This is 
done by the operator n, 

7T = / duiT u , (7) 
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that projects onto colour singlets. Here ft u = e taaQa is a general global U (N c ) 
gauge transformation. The partition function Z is then 



Z = tr (ftp). 



With the colour charge operator 

Q a 



t a 



(8) 



(9) 



the global gauge transformation ft u is in the functional representation: 

^Mv^vM = N u exp [(77* + rj^Q^rii + r] 2 ) + 77^2 - V2V1] ■ ( 10 ) 
Here we have defined 



N u = exp 



1,1/ 1 
-trln-fw + w -1 + 2) 
2 4 



and 



u- 1. 
u + 1 



(11) 



;i2) 



Furthermore, tr = L / tr s trc denotes the trace over all indices, momentum , 

spin and colour, while l s is the unit matrix in spinor space i.e. a 2 x 2-unit 

t a 

matrix (4 x 4 in n = 3). The group element is u = e iaa ~ . Multiplying p and 
u one finds 



(771771 1 up\ rf 2 r} 2 ) = N P N U det(l + fyfi^x 

exp I (77* + 77*) ^ (771 + 772) + ^ - ^771 1 . 

Finally we can then give the expression for the partition function, 
Z = J dutr(ft u p) = J dudet(2)NpN u det(l + ttpQ u ) 



(13) 



J rf M exp|L^tr c ln(l + ^)(l + ^- 1 )|, (14) 



where E = —LN C J p 00 is the vacuum energy. We have here defined 



-0W 



(15) 
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and e* Qj denotes the Nq eigenvalues of the U(Nc) matrix u. On functions of 
eigenvalues the Haar measure du takes the form 



N c 



du = l JJ doti Y[ sin 



2 O-j, ^ 



(16) 



i=l 



This gives 



e -^° / n^Il sin 



2 



X 



exp 



i=l i<j 
N C 



L/"gln(l + ee^)(l + ee" 

^.7=1 



(17) 



In the large A^-limit the integral over U(No) in ( |T7| ) may be calculated by 
the steepest descent method. We take this limit by letting Nc and the spatial 
length L tend to infinity simultaneously keeping their quotient 



Nc 
L ' 



;is) 



constant. We solve this by introducing the density of eigenstates, g, in the 
continuum limit [16, 17] . 



i.e. 



N C r l r a c 

Y^f(a 3 ) - N c / dtf[a(t)} = N c / g(a)daf(a) 

j—l JO J—a c 

We find two distinct phases, depending on the inequality 



(19) 



1 



k Jp e 



> 0. 



(20) 



As long as ( p0|) is satisfied we are in the zero-gap phase, a c = ir, where g 
has support over the whole circle. When (^) is not satisfied we are in the 
one-gap phase. In this case there is a gap in which g lacks support, i.e. an 
interval where the eigenvalues give a vanishing contribution as Nc — > oo. 

In the zero-gap phase In Z is given by a sum over quark — antiquark pairs 



N 



\nZ = -PE - -± \ / In [1 - £(p)£(<z)] . 



(21) 



pjq 



Even though the states are mesonic-like, we do not have confinement since 
InZ oc (Nc/k) 2 = NqL/k, whereas in the confined phase of 't Hooft mesons 
we have (cf. |M) \nZ oc L. In the one-gap phase we find 



-(3E + N 2 \n (sin^) + ^ J ln±(l + £+ ^ + 2£cosa c + 1) 



c 



K Jp Jq 

where we have defined 

x 



In [1 — x(p)x(q)] 



4^ sin 2 ^ 



■1 2 



(22) 



(23) 



Furthermore, the critical angle a c is determined from 



p V V£ 2 + 2£ cos a c + 1 



(24) 



In Figure [l], we plot In Zj (TNq) as a function of the temperature, for k = 
m. The critical temperature is then found as T c « 0.8m. The free energy 
(F = —TlnZ) in the 1-gap phase (solid line) grows as T 2 for large T, and 
is approaching the free energy for free quarks, i.e. without the projection, 
(dotted line), as T — > oo. Whereas the free energy in the 0-gap phase (dashed 
line), continued above the phase transition would grow like T 3 . The high 
temperature behaviour in the 1-gap phase is obtained by an expansion for 
a c 1, yielding a c ~ 4ttk/T. As T — > oo, a c — > 0, i.e. only u — 1 gives a 
contribution, corresponding to no projection. 

In addition to the free energy, it is of interest to compare the entropy, 
and the (renormalized) energy density in the two phases. They are given by 



I A (TlnZ), , 
T 2 d 

T Qf(^Z + (3E ) L , 



(25) 
(26) 



respectively. Their behaviour is similar to the that of the free energy. In 
the interacting case it is well known that a nontrivial chiral condensate ap- 
pears, as shown for m = by Zhitnitsky [[18], 0, and in the general case 
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Figure 1: The free energy for free quarks (dotted line), in the 0-gap phase (dashed 
line), and in the 1-gap phase (solid line), for k = m. All in units of NqT 2 /m. 



by Burkardt f20 |. However, when the interaction only is taken into account 
through the projection, the vacuum is trivial and no chiral condensate can 
appear in this case. Nevertheless, the (renormalized) expectation value of 
ipif), is an order parameter for the phase transition, and thus of interest. We 
find 

(^} = -T^-(\nZ + (3E ) L . (27) 

In the 0-gap phase it is straightforward to evaluate this from (|21|). However, 
in the 1-gap phase we need to define 



a c d 
cot — 1 — — a c = ma, 
2 dm 



where 









I 







lp (^ 2 + 2^cosa c + l) 3 / 2 
is obtained from (p4[), using (0), and (Q). We then find 



T d 

m dm 



x 



x 



This gives 



Ve + 2^cosa c + l 



1 



;i + Oa - (i - 0: 



Nlm{ -~a+ - 



1 



2 k J P 1 + ^ + V^ 2 + 2^ cos a c + 1 
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(28) 



(29) 



(30) 



£ + cos a c 



£ , 2£sin 2 f ^ 
Vf 2 + 2f cosa c + 1 , 



V£ 2 + 2£cosa c + 1_ 
/■ r x( P )x(g) (l + Qa-a-QM 
JpJql- x(p)x(q) V£ 2 + 2£ cos a c + 1 J ' 1 J 

We have been able to show that (^0) is continuous over the phase transition, 
whereas higher derivatives give cumbersome expressions. 

However, in the limit of infinite massive quarks analytical results are 
obtained. As m — > oo, we must also let T c — > oo, such that their quotient 

C = m/T c , (32) 

is kept fixed, in order to keep a finite particle density. Furthermore, we define 
r according to T = T c r. The momentum integrals are now independent of 
p, so we write 

/ - A. (33) 
Jp 

The critical temperature is then determined from 

2 = J P 7^-^ A —l (34) 

In the 0-gap phase we then find 

<« = ^ 2 (l) 2 ^T- < 35 > 

In the 1-gap phase the critical angle is now determined from 

J 1 + 2e^A cos a c + e 2 ^ = ^' + , ] , (36) 
v 1 + k/A 

that after some simplifications give 

*»-*{5 + G)(' + S*^l}- < 37 » 

Notice here the Bose-Einstein distribution in the low temperature 0-gap 
phase, as opposed to the Fermi-Dirac distribution in the 1-gap phase at high 
temperature. It is now a straightforward exercise to show that d 2 / dt 2 {ipip) is 
discontinuous, so that it is a third order phase transition. This is in agreement 
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with the study of Gattringer et al. PTJ concerning static quarks on a line. 



Furthermore, in the corresponding lattice model Gross and Witten [jrjj found 
a third order phase transition, so we conjecture that the phase transition most 
likely is third order also for arbitrary mass. 

ACKNOWLEDGEMENTS 

We are grateful to Ariel Zhitnitsky for proposing the subject, as well as 
stimulating discussions; and to Hans Hansson for interesting remarks, and 
pointing out jTJJ]. J. H.'s research was funded by NFR (the Swedish Natural 
Science Research Council), and D. P.'s by STINT (the Swedish Foundation 
for Cooperation in Research and Higher Education). 



References 



[i] 

[2] 
[3] 
[4] 
[5] 
[6] 
[7 
[8 
[9 
[10 

[11 
[12 

[13 



G. 't Hooft, Nucl. Phys. B B75, 461 (1974). 

C. G. Callan, N. Coote, and D. Gross, Phys. Rev. D13, 1649 (1976). 

T. T. Wu, Phys. Lett. B 71, 142 (1977). 

M. B. Einhorn, Phys. Rev. D14, 3451 (1976). 

R. Brower and W. L. Spence, Phys. Rev. D19, 3024 (1979). 

A. Bassetto and L. Griguolo, Phys. Rev. D53, 7385 (1996). 

A. Bassetto and G. Nardelli, Int. J. Mod. Phys. A12, 1075 (1997). 

A. Bassetto, G. Nardelli, and A. Shuvaev, Nucl. Phys. B495, 451 (1997). 

A. Bassetto, D. Colferai, and G. Nardelli, Nucl. Phys. B501, 227 (1997). 

M. Staudacher and W. Krauth, [hep-th/97091~0l| (1997). 



B. Chibisov and A. R. Zhitnitsky, Phys. Lett. B362, 105 (1995). 

B. S. Skagerstam, Z. Phys. C24, 97 (1984). 

L. D. McLerran and A. Sen, Phys. Rev. D32, 2794 (1985). 



8 



[14] T. H. Hansson and I. Zahed, Phys. Lett. (1993). 

[15] J. Hallin and P. Liljenberg, Phys.Rev. D52, 1150 (1995). 

[16] D. J. Gross and E. Witten, Phys. Rev. D21, 446 (1980). 

[17] G. Mandal, Mod. Phys. Lett. A5, 1147 (1990). 

[18] A. R. Zhitnitsky, Phys. Lett. 165B, 405 (1985). 

[19] A. R. Zhitnitsky, Sov. J. Nucl. Phys. 44, 139 (1986). 

[20] M. Burkardt, Phys. Rev. D53, 933 (1996). 

[21] C. R. Gattringer, L. D. Paniak, and G. W. Semenoff, Annals Phys. 256, 
74 (1997). 



9 



